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ABSTRACT 



We analyze the collisional excitation of the 158 /zm (1900.5 GHz) fine structure 
transition of ionized carbon in terms of line intensities produced by simple cloud 
models. The single C + fine structure transition is a very important coolant of the 
atomic interstellar medium and of photon dominated regions in which carbon is 
partially or completely in ionized form. The [CII] line is widely used as a tracer 
of star formation in the Milky Way and other galaxies. Excitation of the [CII] 
fine structure transition can be via collisions with hydrogen molecules, atoms, 
and electrons. Analysis of [CII] observations is complicated by the fact that it 
is difficult to determine the optical depth of the line. We discuss the excitation 
of the [CII] line, deriving analytic results for several limiting cases and carry out 
numerical solutions using a large velocity gradient model for a more inclusive 
analysis. For antenna temperatures up to 1/3 of the brightness temperature of 
the gas kinetic temperature, the antenna temperature is linearly proportional to 
the column density of C + irrespective of the optical depth of the transition. This 
is appropriately referred to as the effectively optically thin (EOT) approximation. 
We review the critical densities for excitation of the [CII] line by various collision 
partners, briefly analyze C + absorption, and conclude with a discussion of C + 
cooling and how the considerations for line intensities affect the behavior of this 
important coolant of the ISM. 



1. INTRODUCTION 



Velocity-resolved spectroscopy of the 158 /zm fine structure transition of ionized carbon 
([CII]) is currently emerging as a powerful probe of star formation in the interstellar medium 
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of the Milky Way and other galaxies. Due to its relatively low ionization potential, singly 
ionized carbon is the dominant form of the element under a wide variety of conditions. 
[CII] can be observed only from above the Earth's atmosphere due to telluric water 
vapor absorption, except for sources with substantial red shifts. High angular resolution 
(interferometric) observations are not yet available, and velocity-resolved spectroscopy is 
the best avenue for determining what regions are responsible for observed [CII] emission 
features. Such data for Galactic sources have be en extremely rare, and only with the HIFI 



instrument on Hersche 



(IHeyminck et al. 



( Ide Graauw et al. 



2010 ) and the GREAT instrument on SOFIA 



20121 ) have high spectral resolution observations of [CII] become available 



to the astronomical community. 

Unlike the situation for molecular tracers (such as carbon monoxide), for C + emission 
we do not have the possibility of multiple transitions and many isotopologues to allow 
determination of the volume density, temperature, and optical depth, and thus obtain a 
reasonably accurate determination of the column density. In analyzing [CII] data, one 
can compare line profiles with those of CO, its isotopologues, HI, and other tracers, but 
one difficulty in determining the quantity of C + is knowing the optical depth of the single 
transition that is available. In general, the line width of sources is sufficiently large that the 
transitions of 13 C + are blended with those of 12 C + (with the few exceptions discussed in 
§6.11 but which will undoubtedly become more numerous in the near future). The situation 
for absorption is less ambiguous, although in dense foreground clouds the usual assumption 
that all of the C + is in the ground state may not be valid (see §6.3j) . The present analysis 
is intended to show that the line intensities commonly found for 158 /im [CII] emission 
are sufficiently weak that the integrated intensity is proportional to the C + column density 
along the line of sight, irrespective of the optical depth of the line, and hence are in the 
"Effectively Optically Thin" (EOT) limit. 
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2. CLOUD MODEL 

A significant fraction of gaseous carbon is in the form of C + under a wide range of 
astrophysical conditions, including ionized, atomic, and molecular phases of the interstellar 
medium (ISM). Thus, a given line of sight can encompass m ultiple distinct regio ns, as is 



L anger et al. 



20101 ) . Rather 



reflected in the analysis of large-scale [CII] survey data (e.g. 
than consider a highly specific multi-component model from the outset, we consider each 
region as a separate source. Our conclusion, that most of the [CII] emission observed is 
optically thin or in the EOT limit, means that a emission from a given direction can be 
simply decomposed into contributions from the several regions along the line of sight. We 



will thus consider a single cloud, but the results can straightforward 
co mplex situations. T he treatment follows that outl ined by 



by 



Linke et al. 



(119771 ). but which was suggested by 



Penzias 



y be e xtended to more 



Goldreich fc Kwan 



(119751) an d elaborated 



197J) 



We assume that any background energy density in the cloud at the frequency of 
the transition is characterized by a blackbody at temperature T bg . For centimeter- 
and millimeter-wavelength lines, this can be the cosmic microwave background. For 
shorter-wavelength transitions in the immediate vicinity of massive giant molecular clouds, 
this can be the radiation field of the warm dust, which although likely optically thin, can be 
characterized by an equivalent blackbody background temperature that produces the same 
energy density. Note that in general this is not the same temperature as the "background" 
intensity measured in an observation of the source, as discussed in §3.11 and §6.31 

We adopt the escape probability formalism to treat radiative transfer, based on a 
spherical cloud with velocity proportional to radius. However, as the optical depths for 
[CII] are quite modest, the emergent intensity will not be sensitive to the details of the 
velocity field or of the cloud morphology. 



For the present analysis we consider a uniform cloud. Since C + is found in regions 



-5 - 



having temperatures ranging from < 100 K to ~ 10 4 K, this will be an idealization of a 
typical observation. Nevertheless, it does seem that individual spectral features found in the 
observations taken as part of the GOT-C^ survey originate primarily in regions of a sing le 



type 



and physical conditions (ILanger et al. 



2010; 



Velusamy et al. 



2010 



Pineda et al. 



2010l ). We thus adopt a single kinetic temperature T k , C + density n(C + ), and collision 
rate, C. The size of the region is L, so that the C + column density is iV(C + ) = n(C + )L. 

The C + ion has only two fine structure levels in the ground electronic state. The lower 
J = 1/2 level has statistical weight gi = 2. The upper J = 3/2 level has statistical weight 
g u = 4, and lies at equivalent temperature T* = AE /k = 91.25 K above t he ground state. 



19861 ) corresponding 



The measured transition frequency is 1900.537 GHz (jCooksy et al. 
to a transition wavelength of 157.74 /xm. The present discussion is relevant as well to 
13 C + , as long as the hyperfine splitting that is present in this species can be ignored. Th e 



spontaneous decay rate 



Tachiev fc Fischer 



2001 



or the single [C 



Wiese fc Fuhr 



II 15 8 fjm line is 2.3 xlO" 6 s~ l flMendoza 



1983 



20071 ). In £j5]we discuss the critical densities for 



excitation by collisions with different partners. In Table [T] we collect the nonstandard 
symbols that are used in this analysis. 
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Table 1. Nonstandard Symbols Used in Analysis of [CI I] Excitation and Emission 



Symbol 


Explanation or Definition 


rpex 


Excitation Temperature 




Equivalent Temperature (= hu/k) 


rpbg 


Temperature Characterizing Background Radiation Field Internal to Cloud 


rpkin 


Kinetic Temperature 


Ta,so 


Antenna Temperature Produced by Background Source 


G 


Background Term (= [exp(T*/T b9 ) - f]" 1 ) 


K 


Kinetic Temperature Term ( = exp(T* /T hm )) 


X 


Ratio of Downwards Collision Rate to Effective Spontaneous Decay Rate 


9 


Ratio of Statistical Weights ( = g u j 9i) 
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EXCITATION AND EMISSION CALCULATION 



3.1. General two level problem 

The excitation temperature of the transition is defined by the relative populations of 
the upper and lower levels, n u and n h respectively, through the standard relationship 



Qu g— T* /T ex 



ni gi 

Due to the wide range of conditions und er which it is the dominant form of carbon, 



(1) 



collisional excitation of C + by electrons ( 



Wilson fc Bell 1 120021 ) . a tomic hydrogen flLaunav fc Roue: 



by molecular hydrogen (IFlower fc Launay 



Keenan et at. 



1986 



1977 



Flower 



31um Pradhan 



1977: 



1988 



Barinovs et al. 



5 



1991 



20051 ) and 



19901 ) can all be important. 



The upwards and downwards rate coefficients (cm 3 s x ) are related by detailed balance 

Riu/Rui = (g u /gi)e~ T * /Tkm . 



(2) 



For a single collision partner, the collision rates are equal to the rate coefficients times the 
density n of that collision partner, thus 



a 



ui 



R,,,n and C ]u = R, u n 



(3) 



For a region with multiple collision partners, the upwards and downwards rates are the 
sum of the rates produced by each. Within a single region, the kinetic temperatures of 
different collision partners would likely be the same, but this assumption is not necessary 
to calculate the combined collision rates. 



The energy density in the cloud at the frequency of the C + transition is given by 

U = (1 - (5)U(T ex ) + (3U(T b9 ) , 



(4) 



where /3 is the photon escape probability. The background temperature, T b9 , is defined by 
the energy density of the background radiation field at the frequency of the transition of 
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interest. It may be an isotropic radiation field such as the cosmic microwave background 
or the far-infrared background. It might also include a local contribution such as that 
produced by warm dust in a giant molecular cloud with active star formation. For this 
discussion we assume that the C + - emitting region is permeated by this background, 
and that the same background produces the "background" level measured in an observed 
spectrum. The latter latter will not be the case for a cloud that is far from the source of 
the background, which is critical for understanding C + absorption, discussed in §6.31 

For a spherical cloud with a large velocity gradient of the form v oc r, the escape 
probability is given by 



where r is the peak optical depth of the transition. 

Radiative processes include spontaneous emission (rate A u i s -1 ), stimulated emission 
(rate B u iU), and stimulated absorption (rate Bi u U). The stimulated rate coefficients are 
again related by detailed balance through 



(5) 



T 




(6) 



From the relationship between the stimulated and spontaneous downwards rates, 




(7) 



For convenience in dealing with the background, we define 



G 



1 




e T*/T b 9 _ i • 



The rate equation that determines the level populations includes collisional and 
radiative processes, and is 
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n u {A ul + B ul U + C u i) = rn(B lu U + C lu ) , (9) 

With the preceding equations in this section, and substituting into equation [TJ we can 
use equation [9] to obtain an expression for the excitation temperature 



^j.* jrpex _ C u l + 0(1 + G)A u l 

' Gf3A ul + C ul e~ T 'l™ 



(10) 



Defining 



and 



we obtain 



K = e T */ Tkm , (12) 



T*/T ex _ K X + 1 + G (] n) 

6 ~ K X + GK • {V6) 
The quantity $A U \ is an effective spontaneous decay rate, since it is is only those photons 
which escape the cloud that contribute to the downwards transition rate. X is the ratio of 
the downwards collision rate to the effective spontaneous decay rate. Thus, irrespective of 
(3 (and thus of r), for X -> 0, T ex -> T b f , and for X > 1, T ex -> T kin . 

For this two level system, the optical depth can conveniently be described in terms of 
that which would occur if there were no excitation, i.e. T ex = 0: 

r = hB »f C+) . (14) 
ov 

where we approximate the line profile function at line center by 5v~ x , and N(C + ) is the 
total column density of C + . The optical depth in general can be written 

1 _ e -T*/T°* 

T = T °l + (9u/gi)e- T */ T - • (15) 
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We define the ratio of the statistical weights by 




(16) 



and using equation [13] the optical depth can be expressed as 



X(K -1) + K 



(17) 



T = T 



X(K + g) + K[l + G(l+g)} 



3.2. Antenna Temperature and Specific Intensity 



If, as discussed in ^J2j the background radiation density in the cloud is characterized by 
a temperature T b9 that is also responsible for the "background" temperature in a spectral 
line observation, the antenna temperature (above the background) for a uniform source 
that fills the entire antenna beam pattern is given by 



This will not be the case for clouds for which the background radiation source subtends only 
a small solid angle, as for clouds arbitrarily located along the line of sight to the source. 
This situation needs to be treated somewhat differently, as discussed in §6.31 

The antenna temperature is proportional to the power collected by the telescope per 
unit bandwidth, which is proportional to the collecting area of the telescope (the effective 
area A e in antenna terminology) and to the normalized response convolved with the 
brightness distribution on the sky. For a uniform source that fills the entire antenna beam, 
the integral becomes the product of the antenna solid angle Q a and the source brightness. 
For a diffraction limited antenna, A e Q a = A 2 . The result is that the antenna temperature is 
proportional to the specific intensity I v (W m~ 2 sr _1 Hz -1 ), through the relationship 



i* 



AT A = ( 



I e T*/T b s _ 



T )[l-e"1 



(18) 



AT 




(19) 
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where AI U is the specific intensity of the line above that of the background. For the 
[CII] line, the numerical value of the conversion constant in equation [TH] is 9xl0 14 for 
the specific intensity in the units given above. Many observations of the [CII] as well as 
other sub millimeter spectral lines are carried out by low spectral resolution systems which 
measure only the integrated intensity. For conditions as described above, there is a simple 
relationship between the intensity and the antenna temperature integrated over the spectral 
line in question. For [CII] the relationship in generally-used units is 



With equations [T21 [T71 and [TBI we can calculate the observed antenna temperature for 
arbitrary background and kinetic temperature as well as excitation rate and C + column 
density. We now consider some limiting cases that are also of practical interest. 




(20) 



3.3. Limiting Cases 



3.3.1. Optically Thin Emission 



For low opacity (r<l and ~ 1), we have 



X 




(21) 



and find that 




(22) 



Substituting equation [T7] for the optical depth, we obtain 



AT A (thin) = T 



X[1-G{K-1)] 



TO , 



(23) 



X(g + K) + K[l + G(l + g)} 
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in which the dependence on the total C + column density and the excitation are clearly 
separated. The behavior of the antenna temperature as a function of collision rate (through 
X) for three kinetic temperatures, with and without background, is shown in Figure [TJ 

20 



1.5 



1.0 

o 
< 

M 
O 

— ' 0.5 



0.0 



-0.5 

-2.0 -1.0 0.0 1.0 2.0 

log(X = C ul /A ul ) 

Fig. 1. — Antenna temperature for optically thin C + emission for three kinetic temperatures. 
The antenna temperature is normalized by the no-excitation optical depth To (equation [T7|) . 
The results with no background radiation are indicated by solid lines, while those for a 
background temperature of 40 K are indicated by dashed lines. This very large value of the 
background indicates clearly how the observed antenna temperature (above the background) 
is reduced by the presence of the background. 




In many situations, the background is negligible, and we can simplify the preceding 



13 



equation to yield 



ATA(thin;no background) = T*—^ — — ^ — , (24) 

which can also be expressed as 

AT A{thi n; no 6„c W ) = * £ A - "f > . (25) 

zitkv 1 1 + £-(1 + dv 

Substituting the relevant quantities for the [CII] 158 urn. line, we obtain 

2 4x1 0~ 6 N(C + ) 

AT A (thin; no background) = 3.43 x 10- 16 [l + 0.5e 9L25/TfcI "(l + - — , (26) 

C u i ov 



where the line width is in km s . ICrawford et al. 



(]1985[ ) give essentially equivalent 



expressions for the specific intensity for the [C + ] line, but not considering any background 
radiation. 



3.3.2. Optically Thin Subthermal Emission 

For optically thin emission ((3 ~ 1), subthermal means that C u i/A u i < 1. Equation 1221 
in this limit becomes 

AT A {thin; subthermal) = — -^[1 - G(K - 1)] . (27) 

K A ul 



In this limit, we also find 



t( subthermal) = , (28) 

v ' l + G(l + g) v ' 



and we can write 



T* X 1 G(K 1) 

ATA(thin; subthermal) = — - — — - — — -r . (29) 

This result can also be derived from equation [23] in the low excitation limit, and is 
convenient, as Tq is independent of excitation and is proportional to the total column density 
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of C + , as seen from equation [TH Equation [29] reflects the fact that in the low-excitation 
limit, the population distribution is determined by the background. The fraction of the total 
population in the ground state is (1 + G(l + <?)) _1 , and the observed antenna temperature 
is proportional to the rate of excitation from the ground state. 

Substitution of definitions of the radiative and collisional rates allows equation [2H] to 
be written as 

he 3 1 — G(K — 1) 
AT A (thin;subthermal) = - C lu - ^ ' N(C + ) , (30) 

cmkv l ov 1 + G(l + g) 

where Sv is the line width in cm s _1 . The value of the deexcitation rate coefficient will 

depend on the collision partner and may itself be dependent on the kinetic temperature. 



3.3.3. Optically Thin Thermalized Emission 

At high densities, C u i 3> A u[ and in the corresponding limit X > 1 we find 

AT A (thin;thermalized) —T* ^~ T o ■ (31) 

For thermalized optically thin emission, the background reduces the observed antenna 
temperature by a factor 1 — G(K — 1), which is closer to unity than the analogous reduction 
factor for subthermal excitation. 



3.3.4- Optically Thick Emission 



For high opacity (t»1 and ~ r 1 1), we have 



X = , (32) 

A u i 



and find that 

< c ^)\l-G(K -1)] 



AT A ( thick) = T* y . (33) 
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Equations [22] and [33] are very similar but not identical, as the optical depth appears in the 
denominator for optically thick emission. 

The behavior of equation [33] is illustrated in Figure [2] which shows the antenna 
temperature for the C + line as a function of X = C u i/(3A u i for three different kinetic 
temperatures, assuming optically thick emission. The background is taken to be zero. The 
antenna temperature is approximately a linear function of X, for X as large as unity, or C u \ 
= (3A ul . 

As seen in equation [33] the effect of the background is to reduce the intensity by 
a "background factor" equal to 1 — G(K — 1), independent of the value of X. This is 
illustrated in figure [3] which shows even for T hg = 60 K , the weakening is a modest factor 
of about 2, while for T bg = 90 K, it is a very significant factor of 7. 



3.3.5. Optically Thick Subthermal Emission 

Subthermal excitation for optically thick emission with r>l means (from equation 
[TBI that X = C u i/(3A ul = C u it/A u1 <C 1, in which limit equation [331 becomes 

AT A (thick; subthermal) = ^%^[1 - G(K - 1)] , (34) 

K Aui 

which we see is identical to equation [27] Substituting the expression for tq (in the limit X 
1), we obtain an expression identical to equation [29] This confirms that in the limit of 
weak excitation, J < 1, the observed intensity is linearly proportional to the total column 
density even though the optical depth is large. 
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0.0 
-1.0 

-2.0 -1.0 0.0 1.0 2.0 

Log (X = C ul //SAJ 

Fig. 2. — Antenna temperature of C + 158/im transition in the optically thick limit as a 
function of X = C u i/j3 u i, for three kinetic temperatures. In this limit, (3 = r" 1 , so that 
the antenna temperature is linearly proportional to the optical depth and thus the column 
density as long as X < 1, despite the fact that r» 1. There is no background radiation in 
this example. 

4. General Behavior of the Antenna Temperature 

The preceding discussion has indicated the behavior of the C + fine structure line in a 
number of limiting cases. Unfortunately, the dependence of the excitation on the optical 
depth and the fact that r is itself a function of the degree of excitation are significant 
obstacles to a general analytical solution. We have thus resorted to numerical solution of 
the rate equations including the effect of radiative trapping and collisional excitation and 
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2.0 



^ 1.0 
< 

0.0 
-1.0 

-2.0 -1.0 0.0 1.0 2.0 

Log (X = Cy/BAJ 

Fig. 3. — Effect of background on the antenna temperature of C + 158 /im transition for 
optically thick emission. The background is characterized by a blackbody temperature T& s , 
while the kinetic temperature is 100 K. 

deexcitation as described in §3.11 This allows us to confirm the limits on the conditions 
under which the antenna temperature is linearly dependent on the C + column density. 

We first consider the behavior of ATa as a function of r . As given by equation the 
zero-excitation optical depth is proportional to the C + column density per unit line width. 
For the [CII] line, this is given by 

r = 7.49 x 10- ls N(C + )/5v(km/s) , (35) 

where Sv is the line width. Figure H] shows AT4 as a function of r for a representative 
situation in which there is no background radiation. 
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log (t ) 

Fig. 4. — Variation of the antenna temperature of the C + 158 fim transition as a function 
of the zero-excitation optical depth To, which is proportional to the column density of C + 
as given by equation |35j The kinetic temperature is 100 K, and the different curves are for 
different ratios of C u \j A u \. 

As one approaches thermalization, the antenna temperature becomes independent of 
the collision rate, since the level populations approach their values characteristic of the 
kinetic temperature. The antenna temperature is also limited by the brightness temperature 
equivalent of the kinetic temperature, which is 61.2 K for the kinetic temperature of 100 K. 

Over much of the range of tq the antenna temperature is a linear function of To. For 
subthermal excitation, is proportional to r and to the collision rate, and this situation 
can apply even to cases in which r 3> 1, as indicated by the preceding discussion. 
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From equation [T71 for no background 

X(K-1) + K 
T = T0 X(K + g) + K ' (36) 

The actual optical depth r thus drops from To in the limit X <C 1 to tq{K — 1)/(K + g) in 

the thermalized limit with X ^> 1. For C + at a kinetic temperature of 100 K, the reduction 

is a factor ~ 3, so that the actual optical depth is always within this factor of r . It is only 

as the population of the upper level becomes significant as a result of collisions (possibly 

enhanced by radiative trapping) that the optical depth starts to differ from r , which also 

defines the breakdown of the linear relationship between the antenna temperature and the 

column density. 

It is apparent that there are two limits for this linear behavior. For thermalized 
emission, the excitation is fixed and the optical depth is a constant times the column 
density. For r = 0.5, the antenna temperature is ~ 0.8 times the extrapolated value for 
optically thin emission, tq = 1.5 might thus be considered a limit to the linear regime, 
corresponding to an antenna temperature ATa = 24 K (for T kin = 100 K). 

For non-thermalized emission, the limit is on the quantity C u iTq/A u i rather than just 
on the optical depth. For a fixed value of C u iTq/A u i, the deviations from linear behavior 
increase as To increases. For C u i/A u i < 1, a limit C u it /A u i < 1 ensures no more than 40% 
deviation, and more stringent limit C u it /A u i < 0.5 guarantees linear behavior to within 
20%. 

In Figure |5] we show the dependence of the antenna temperature on the actual optical 
depth, which is itself computed from the solution of the rate equations. While very similar 
to Figure HI there are some subtle differences due to the variation of r resulting from the 
changing excitation temperature. As discussed above, the transformation from r to r 
depends on C u ij A u \ and on 0. This is seen in the deviation of the horizontal positions of 
the triangles, which are defined by fixed values of N(C + ) and thus r . As the excitation 
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increases, the optical depth drops for a given C + column density. 



2.0 I i i i i | i i i i | i i i i | i i i i | i i i i | i 




log (t) 



Fig. 5. — Variation of the antenna temperature of the C + 158 /im transition as a function of 
optical depth r. The kinetic temperature is 100 K, and the different curves are for different 
values of C ul /A ul . 

In the subthermal limit r = To, so the antenna temperature is proportional to both of 
these quantities. The demarcation of the limit in which AT a is proportional to r is very 
similar to that defined by r . For thermalized optically thin emission, r is a factor of a few 
lower than r for reasonable values of the kinetic temperature. In consequence, the limit of 
the linearly proportional relationship is very similar to that derived for ro. Following the 
above discussion, we find that for C u \j A u \ > 10 and r < 0.5 we have a linear relationship 
to within 20%. For subthermal excitation with C u \r j ' A u \ < 1, the combined limit C U \T j A u i 
< 1 ensures linear behavior to within 40%. The upper limit on the observed antenna 
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temperature (for T km = 100 K) is again ~ 20 K, which more generally corresponds to ~ 
1/3 the difference between the brightness temperature of the kinetic temperature minus 
that of the background temperature. The nonlinear regime starts at slightly lower values of 
t for lower kinetic temperatures and for significant background. A more conservative and 
broadly applicable limit for sub thermal and thermalized emission is given by C U \T / ' A u \ < 
0.5 is appropriate to broadly ensure proportionality between column density and antenna 
temperature to within 20%. This can be considered to be the limit of the effectively 
optically thin (EOT) approximation. 



5. COLLISION RATE COEFFICIENTS AND CRITICAL DENSITIES 

Excitation of the [CII] transition can be through collisions with electrons, hydrogen 
atoms, and hydrogen molecules, since C + can be a significant form of carbon in a wide range 
of conditions in the interstellar medium. For each type of collision partner, the deexcitation 
rate C u \ is the deexcitation rate coefficient R u i multiplied by the collision partner density 
n. The critical density n cr is defined by the collisional deexcitation rate being equal to the 
effective spontaneous decay rate. For optically thin emission, the latter is just A u i, while 
more generally it is given by j3A u i. We thus have 

(3A ul 

n cr = — — . (37) 

tiul 

From equation [TTJ, this leads to the definition 

X = — . (38) 

Tt cr 

For optically thin emission with no background radiation, we see from equation [T3l that 

T"M = jJ* 9 • (39) 

For C + emission from an atomic hydrogen cloud or warm PDR region, T* ~ T kin , so that 
T ex {n cr ) ~ 50 K. 
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Collisions with elec t rons have been analyzed by many authors ( IKeenan et al. 111986 



Blum fe Pradhan 



1991 



1992 



Wilson fe Bell 



20021 ). Over a wide range of temperatures 



from ~ 100 K to 20,000 K, the deexcitation collision rate coefficient as a function of the 
electron temperature can be written 



R,Ae 



8.7 x 10~ 8 (T e /2000) 



"°- 37 cmV 1 



(40) 



where T e is the electron temperature. 



Collisions between C + and hydrogen ato ms have been st u died 



( 119771 ) and more recently in great detail by 



Barinovs et al. 



by 



Launay fe Roueff 



fl2005f ). The collisional 



deexcitation rate coefficient as a function of temperature from the most recent analysis is 



R u i{H°) = 4.0 x 10- n (16 + 0.35T - 5 + 48T" 1 ) 



s — 1 

cm s 



(41) 



A satisfactory fit over a the range 20 K < T km < 2000 K is given by 



RuiiH") = 7.6 x 10- 10 (T fein /100) ai4 cmV 1 



(42) 



The situation for collisions with H2 molecules is more complicated due to the 
existence of two different spin modifications of the molecule, the relatively uncertain ortho- 
to-para-ratio, and the rotational s tates of the molecule that may be populated. The earlies t 



calc ulations by 



and 



Flower 



Ch 



u fe Dalgarno 



( 119751 ) were impr oved upon by 



Flower 



'lower fe Launav 



(119771) 



(119901 ) . From his cooling rates 



(119881 ) . The situation was reviewed by 
for C + excited by collisions with H 2 compared to excitation by collisions with H°, we infer 
that over a broad range of temperatures, the collision rate coefficients for collisions with 
molecular hydrogen are approximately half those for collisions with atomic hydrogen. There 
are minor differences for the two spin modifications, but the main reason for the smaller 



rate coefficient for H2 is its lower velocity at a given temperature, as the cross sections for 
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H2 - C + collisions are very similar to those for H° - C + collisions. Over the temperature 
range appropriate for molecular clouds we thus adopt 

R u i{H 2 ) = 3.8 x l(T 10 (T fci 7i00) a14 cmV 1 . (43) 

Table [2] gives the resulting critical densities for different collision partners. We see that the 
critical densities for electron excitation are about a factor of ~ 50 to ~ 500 smaller than for 
excitation by atomic or molecular hydrogen. 

The expression n/n cr can be substituted for X in any of the equations in §3.31 In 
general, for densities below the critical density (X < 1) the excitation is sub thermal, and 
for n ^> n cr (X ^ 1) we achieve thermalization. In the event that a single region has 
multiple collision partners (e.g. H 2 and H°), the value of X is determined by the sum 
of the downwards collision rates compared to the spontaneous decay rate. In practice, 
while the C + emission from a given line of sight may include contributions from primarily 
molecular, atomic, or ionized regions, these are likely to be physically distinct, and have 
quite different kinetic temperatures. If the combined emission is optically thin, the total 
is the sum of that from each region. The total emission can still be computed using an 
expanded version of equation [23] with the C + column density and line width, together with 
the kinetic temperature and background for each region. 
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Table 2. Critical Densities for [CII] 158 /iin Fine Structure Line (cm 3 ) 



Temperature 




Collision Partner 




(K) 


e 


H° 


H 2 


20 


5 


3800 


7600 


100 


9 


3000 


6100 


500 


16 


2400 


4800 


1000 


20 


2200 


4400 


8000 


44 


1600 


3300 
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6. Application to Observations 



6.1. Galactic Emission 



New observational facilities are rapidly increasing the quality and quantity of [CI I] 
observations to the point where making a complete review of [CII] fine structure emission is 
not feasible We here focus on a very brief overview indicating some of the cases where the 
previous discussion can be applied and some where conditions make it unlikely to be helpful. 
An obvious concern is that the very simplified model, particularly the assumption of single 
kinetic temperature and collision rate, means that realistic lines of sight will require more 
complex, multicomponent models to model accurately the observed emission. Nontheless, 
there appear to be many situations in which general considerations outlined above can be 
applied. 

There have been a number of large-s cale surveys of [CHI in the Galaxy. The 



FIRAS instrument on the COBE satellite (IBoggess et al. 



19921 ) made essentially all-sky 



spectrally-u nresolved measurements of the [C II] line with 7° angular resolution. As 



reported by 



Fixsen. Bennett. &: Mather 



(119991 ). the intensity varied from a maximum of ~ 
150 K kms -1 in the inner galaxy, to ~ 40 K kms -1 in the outer galaxy. With reasonable line 
widths, the peak equivalent antenna temperature is < 10 K. Thus, the effectively optically 
thin (EOT) approximation should be applicable, subject to the significant uncertainty and 
variation in the range of conditions and possible severe beam dilution of small sources in 
the very large beam. 



Nakagawa et al. 



( 119981 ) used a balloon-borne instrument having a 20 cm effective 
telescope diameter producing a beam size of 12'. 4 (FHWM) and a Fabry-Perot spectrometer 
yielding a velocity resolution of 175 km s _1 . The survey covered a large fraction of the 
inner Milky Way but without spectrally resolving the C + emission. Typical strong peaks 
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given in their Table 2 have intensities I ~ 10~ 3 erg s~ 1 cm~ 2 sr~ 1 . From equation 1201 this 
is equivalent to a peak antenna temperature of 14 K assuming a line width of 10 km s _1 . 
This suggests that the effectively optically thin approximation is satisfactory and that the 
observed intensity is proportional to the column density of C + . 

The availability of sub-km s _1 r esolution with the he terodyne instrument HIFI 



(Ide Graauw et al. 



20101 ) on Herschel flPilbratt et al 



20101 ) . together with the 12" angular 



resolution afforded by its 3.5 m diameter telescope, enabled the GOT-C^ Open Time Key 
Project. This survey of the Galactic Plane sampled some 900 lines of sight covering 360° 
in longitude, but not targeting specific region s and, of course, being highly undersamp led 



in angle. The limited data released to date (ILanger et al. 



2010 



Pineda et al. 



2010; 



Velusamv et al. 



2010. 



20121 ) reveal complex spectra with multiple kinematic components. 
Integrated intensities are less than 10 K kms -1 , and peak antenna temperatures rarely 
exceed 3 K. The spectra from the GOT-G^ project shown in figure [6] illustrate the typical 
weakness of the [CII] lines observed in the Galactic Plane when no effort has been made to 
select regions with ongoing star formation. Figure shows the distribution of peak main 
beam brightness temperatures from the GOT-C^ survey; the most commonly-found values 
of T m b, within a few tenths of a K of zero, are due to the radiometer noise. The low intensity 
of the [CII] emission suggests that the EOT approximation can be used to facilitate data 
analysis. 

The situation is quite different if angularly-resolved bright PDR regions associated 
with massive young stars are observed. This was evident from observations made with 
incoherent detectors that while not fully resolving the line profiles indicated very strong 



emission and for reasonable val ues of line width, 



include observations of Ml 7 by iMatsuhara et al. 



a rge a ntenna temperatures. Examples 



(1l989f ) with a 3'. 7 beam, who found 



a peak intensity of 340 K kms 1 . For a line width of 10 km s 1 this corresponds to 
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Fig. 6. — Selection of [CII] spectra from the GOT-C^ survey. These spectra were taken at 
quasi-uniformly spaced directions in the Galactic Plane. The vertical axis is the antenna 
temperature corrected for the Herschel main beam efficiency 



an antenna temperature of 34 K, which is nearing the limit for the effectively optically 
thin approximation. If the kinetic temperature in the PDR region is in excess of 150 
K, as predicted by some models, the approximation might still be satis factory. Recent 



observations of resolve d PDR regions include S140 (IDedes et al. 



HIFI and Cepheus B (IMookerjea et al. 



2010) observed with 



20121 ) observed by the GREAT instrument on 



SOFIA. These both yield moderately strong lines, which are close to the limit of the EOT 
approximation. 
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100000 




T mb /K 



Fig. 7. — Histogram of main beam temperatures from the Herschel GOT-C^ survey. The 
obvious peak, centered around zero K, is due to the radiometric noise from the HIFI het- 
erodyne instrument. The distribution of the noise is quite well represented by the Gaussian 
shown (which was fitted for negative values of the main beam temperature). The [CII] emis- 
sion, indicated by the shaded region, is the excess above this noise, and is evident for T m b 
between ~ 0.5 K and 8 K. Very few spaxels are present with T m {, > 8 K. 



Observations of Orion by I Crawford et al. 



(I1986I ) yielded a peak intensity of 5.5x10 



erg s _1 cm _2 sr _1 , corresponding to a peak antenna temperature of 150 K if the line width 
is 5 km s -1 (as revealed by subsequent heterodyne observations). These authors concluded 
that the [CII] emission is optically thick wi th r ~ 1. The 10 km s _1 velocity resolution of 



the Fabry-Perot spectrometer employed by 



Stacey et al. 



(1991) did not quite resolve the 
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[CII] line. The peak intensity measured is 4.1xl0 -3 erg s _1 cm _2 sr _1 , which for 5 km s _1 
line width corresponds to an antenna temperature of 120 K. The velocity resolution was 
sufficient to resolve the F = 1—^0 component of the [ 13 CII] fine structure line and these 
authors conclude that while the [CII] optical depth at the emission peak may approach 2, 
the emission throughout most of the region mapped is optically thin. 

The Orion PDR region is considered to have densities above 10 5 cm -3 , at least an 
order of magnitude larger than the critical density for collisions with H° or H2 (Table |2]). 
With C u i/A u i > 100, the transition will be thermalized. The background continuum is 
equivalent to an antenna temperature of < 10 K and has only minor effect on the emergent 
line intensity. The two analyses discussed above are consistent in suggesting that the optical 
depth is ~ 1 in this PDR. Since we are not sure to be optically thin, the discussion of §1] 
is relevant. The optical depth can be determined if the kinetic temperature of the emitting 
region can be determined (e.g. from model of PDR thermal balance) or inferred from the 
[ 13 CII] line if the isotope abundance ratio is known. 



he first heterodyne observations of the [CII] line were made by 



Boreiko. Betz. fc Zmuidzinas 



( 119881 ) who measured a line width of 5 km s 1 and a peak antenna temperature of 86 K 



with a 55" (FWHM) beam. While somewhat lower than that inferred from incoherent 



measurements, the ability to resolve the F = 2— >1 hyperfine component o 



indicated an optical depth r([ CII]) ~ 5 



Boreiko. Betz. fc Zmuidzinas 



the [ 13 CII] line 



( 11990( 1 observed 



a number of bright PDR sources with 0.8 km s _1 velocity resolution and 55" (FWHM) 
beam size, finding peak antenna temperatures between ~ 30 K and ~ 70 K. If the kinetic 
temperatures in the PDRs are in excess of 100 K, their high densities and resulting 
thermalized emission suggests that the optical depth may be modest, but this is not certain. 
This study also suggested that a significant portion of the [CII] emission may arise in 
ionized gas, making it more difficult to determine the C + column density. 
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The Herschel Guaranteed Ti me Project WADI 



of the HII region NGC 3603 MM2 flOssenkopf et al 



report ed velocity-resolved observations 



20111 ). The 12 C + fine structure line has 



a peak antenna temperature of ~ 20 K. Based on observations of the F = 1—^0 component 
of the 13 C + transition, the optical depths may be be small or moderate, depending on the 
13 C + / 12 C + abundance ratio. 



Graf et al. 



(120 121 ) show a striking 



The GREAT observations of NGC2024 by 
self-reversal in the 12 C + line. This strongly suggests that there is low excitation, optically 
thick gas in front of the bright PDR emission, and that the EOT approximation as 
described above cannot be directly applied. The assertion that the relative intensities of 
the 13 C + hyperfine components cannot be fit by any model in which saturation modifies the 
observed ratio relative to the optically t hin ratio may be due to a problem with the assumed 



relativ e line st rengths (lOssenkop 



(120 121 ) and of 



Ossenkopf et al. 



20121) . Despite this complication, the data of 

IT 



Graf et al. 



(120121 ) dramatically illustrate the enormous potential value 
of high resolution spectroscopy in determining the optical depth of the [ 12 CII] line and 
understanding its relationship to luminous PDRs in the Milky Way as well as in external 
galaxies. 



6.2. Extragalactic Emission 



[CII] emission has been recognized as an important tracer of star formation 
throughout the universe, to a significant degree because of its great luminosity, which 
i s typically between . 1% to 1% of the far-infrared luminosity o f star-forming galaxies 



(Madden et al. 



1993 



Malhotra et al. 



2001 



Stacey et al. 



20101 ). The distribution of 



[CI II within external galaxies ha s been studied to a limited extent (jPoglitsch et al. 



1995 



Rodrigez-Fernandez et al. 



20061 ). but the capability o 



has allowed much improved studies of its distribution (IMookerjea et al 



the Herschel PACS instrument 



20111 ). in part to 
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determine the fraction of the line emission originating in molecular, atomic, and ionized 
regions. This ubiquity of C + with the varying contributions of these different components 
of the interstellar medium to the total energy radiated by C + i n different galaxies m akes 



calibration of the [CI 



de Looze et al. 



versus star-formation rate (discussed by 



Boselli et al. 



fl20D2h and 



(1201 ll )) a difficult task. 



From the present discussion we can certainly see that emission from dense PDRs 
associated with massive young stars is likely thermalized and possibly optically thick. The 
emission thus varies as n°, and as N a , with < a < 1, in contrast to the subthermal, 
optically thin emission from the diffuse atomic component of the interstellar medium, which 
varies as n 1 and N 1 . The discussion of Galactic emission makes it reasonable that a fraction 
of the [CII] emission will be optically thick, but this will likely be a small fraction of the 
area, if not the total intensity, in most sources. The problem of the line opacity is distinct 
from that of possible absor ption by dust, which has been suggested as bei ng significant at 



the wavelength of [CII] by 



Papadopoulous. Isaak. fc Van der Werf 



f|2010h . The fact that 



there is as good a correlation as is observed between tracers of star formation and [CII] 
emission is likely a reflection of the fact that star formation is associated with PDRs, HII 
regions, and relatively massive atomic clouds. Significant variations would not be very 
surprising, and the situation may only become more complex when studied with higher 
angular resolution data as will soon be available from ALMA for red-shifted galaxies at z 
> 1. 



6.3. Absorption by Foreground Clouds 

The clouds modeled in §3] are assumed to be permeated at the frequency of the 
spectral line of interest by a blackbody radiation field characterized by T hg . As a result, the 
excitation temperature of the two level system is limited to the range T bg < T ex < T km , and 
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there can be no absorption unless T km is less than T bg , which is not physically plausible. 
The situation is different if we have regions with disparate background radiation fields. This 
can occur if we have a foreground cloud (with no particular background radiation field) 
along the line of sight to a strong source of continuum and line emission associated with 
massive star formation in a GMC In this case, the excitation temperature of [CII] in the 
foreground cloud can be as low as the temperature of the CMB , and it is certainly plausible 



that we see this cloud in absorption. 



Falgarone et al. 



(120 lOl ) observed CII absorption 



(line intensity less than that of continuum) towards the Galactic HII region DR21. The 
distinction between the background and foreground clouds is facilitated by the line of sight 
velocity difference resulting from Galactic rotation and the long path lengths involved. 

The Herschel Guaranteed Time Key Project (GTKP) PRISMAS (M. Gerin, P.I.) 
targets strong continuum sources specifically to probe absorption of many species along the 
line of sight. The GTKP Project HEXOS (E. Bergin, P.I.) has detected [CII] absorption 
towards Sgr B2. We denote the antenna temperature (referred to a single sideband) 
produced by the background source Ta, S o- Using equation [181 we find that to get absorption, 
T ex < T*/ln(l+T*/T A , so ). 

i — ii — in 

Ta,so varies from ~ 2 K to ~ 10 K (IGerin ll2012al Jbl). The maximum C + excitation 
temperature in the foreground cloud to get absorption ranges from ~ 25 K for Ta, S o = 
2.5 K to ~ 40 K for T A , so = 10 K. For the sources with stronger background antenna 
temperatures, the C + line is generally seen in absorption, since the excitation temperature 
in the line of sight clouds can well be less than 40 K. 

Equation [13] allows us to obtain an upper limit on the density in the foreground 
cloud that is seen in absorption. Assuming that the background radiation field in the 
foreground cloud is zero, we find T ex = T*\T* /T km + ln(l + 1/X)]" 1 . Taking the kinetic 
temperature of the foreground cloud to be 100 K and assuming that the [CII] line is 
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optically thin and excited by collisions with atomic hydrogen, we find from $51 that n(H°) < 
200 cm" 3 for T AySO = 2.5 K and n(H°) < 1000 cm' 3 for T AySO = 10 K. These densities are 
generally consistent with what is known about these relatively diffuse foreground clouds. 
The non-zero excitation of the C + means that determining the column density from the 
calculated spectrum is not as straightforward as for absorption lines from species for which 
all of the population can be assumed to be in the ground state. Equation [35] can be used to 
calculate the optical depth if the excitation temperature or density is known, and equation 
1351 can then be used to derive the C + column density in the foreground cloud. 

If the C + line is assumed to be optically thick, we can use the minimum temperature 
to obtain T ex and then from equation [13] obtain X. An estimate of the optical depth is 
required to derive the collision rate and thus the density, but this can yield both the density 
and C + column density of the foreground cloud. 



7. C+ Cooling 



The cooling rate in erg s 1 cm 3 from the single C + fine structure transition can be 



written ([Goldsmith fc Langer 



19781 ) 



A = (A ul n u + B ul n u U{T h9 ) - B lu mU{T ha ))(5h 



v 



(44) 



where the upper and lower level C + densities, n u and ni, that result from the combination 
of collisional and radiative processes, are in cm -3 . The first term reflects the spontaneous 
emission that escapes the cloud, the second term is the radiation produced by stimulated 
emission by the background radiation field (at temperature T bg ) that escapes from the 
cloud, and the third term is the radiation absorbed from the background radiation field. 
The cooling rate can be expressed as 

e (T*/T°*) _ l 



A = n u A ul f3hv[l 



e (T*/T b a) _ 



(45) 
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A = Ci/ ^ 777 777 -n(C + )gA ul hu . (47) 



with the term in square brackets being the background correction term to the usual 
expression for the cooling rate. Using the terminology from §3.11 we can write this as 

A = x^7W7W7W] 9flAMC+)h " • m 

where n(C + ) is the total C + density. For optically thin emission, X = C u i/A u i and (3 = 1. 
We see that in this limit equation H6] becomes 

<fc[l-G(K-l)] 
^.( g + K)+K[l + G(l + g)} 

This form is identical to that of equation [221 except with gA u [hv replacing T*r . This 
reflects the fact that for a system with a single transition, the antenna temperature (as well 
as of course the specific intensity) and the cooling are proportional. Thus, the previous 
considerations of thermalization and limits on proportionality between C + column and 
emission intensity are all applicable. 

Equation H3 for the case of no background (G = 0) can be written 

A = A ul hvn(C+)[l + -(1 + ^)]" 1 , (48) 

which is, of course, very similar to equation |25j As an example, we may consider a PDR in 
which the gas temperature is on the order of 100 K, and the hydrogen is largely molecular. 
If the collision rate is much less than the spontaneous decay rate, equivalent to the density 
being well below the critical density (~ 6000 cm -3 ; Table [2]), the cooling rate is proportional 
to the C + density times the H 2 density. Adopting the rate for collisions with H 2 discussed 
in we can write the cooling rate as 

A = 9.6 x 10- 24 e- 91 - 25 / Tfcl "(^^) 014 n( J ff 2 )n(C' + ) er^W 3 . (49) 



This expression is the form for optically thin subtherm al cooling that is almost always 
given in treatments of the neutral ISM (e.g. iDraine 1 1201 ll ). More generally, however, 
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especially considering [CII] line cooling of bright, dense PDRs, we must consider the 
situation in which the density is high enough that collisional deexcitation is significant, and 
that the transition may even be thermalized. As seen from equation H3, the cooling per C + 
ion is reduced in this situation. A similar effect is produced by finite opacity. As discussed 
above, the general case has to be calculated from solution of the nonlinear rate equation, 
but the results can be obtained by scaling those for the antenna temperature presented in 

SI 

8. CONCLUSIONS 

We have analyzed the radiative transfer and collisional excitation of the [CII] 158 fim 
fine structure transition. Using an idealized cloud model based on the escape probability 
approach, we have presented general expressions for the [CII] antenna temperature that 
are straightforward to convert to the specific intensity of the fine structure line. Much 
of this analysis is applicable to other two level systems, and in may cases to the ground 
state transition of multilevel systems (such as, for example, the water molecule). In various 
relevant limits, analytic expressions can be obtained, which are confirmed by numerical 
solution of the rate equations including trapping appropriate for spherical cloud large 
velocity gradient model with v oar. In the weak excitation limit defined by X = C u i/(3A u i 
<C 1, the Effectively Optically Thin (EOT) approximation indicates that since every photon 
produced by spontaneous emission ultimately escapes the cloud, the antenna temperature 
is proportional to the C + column density even if the optical depth r>l. 

The strongest lines for which the EOT approximation is valid have ATa < 0.3 times 
the brightness temperature of the kinetic temperature minus that of the background 
temperature, or typically 20 K to 30 K for [CII] line emission from clouds in the ISM. The 
[CII] fine structure line emission from well-resolved PDR regions may be too strong to be 
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in this limit, but interpretation of much of the [CII] emission traced by large-scale surveys 
of the Milky Way, and which is likely dominant in a large fraction of the area of external 
galaxies, will be facilitated by the EOT approximation. [CII] cooling from the single 
fine structure is proportional to the antenna temperature, and consideration of collisional 
deexcitation and finite optical depth will reduce the cooling per C + ion relative to that 
obtained in the case of subthermal and optically thin emission. 
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